arXiv: 1506.08645vl [astro-ph.HE] 29Jun2015 


Energy bursts from deconfinement in high-mass twin stars 


D. E. Alvarez-Castillo^’^’“, M. Bejger^, D. Blaschke"’^, P. Haensel^, L. Zdunik^ 
°‘Institut Fizyki Teoretycznej, Uniwersytet Wroclawski, pi. Maxa Borna 9, 50-204 Wroclaw, Poland 
^N. Copernicus Astronomical Center, Polish Aeademy of Sciences, Bartycka 18, PL-00-716 Warszawa, 

Poland 


ABSTRACT 

We estimate the energy reservoir available in the deconfinement phase transition induced 
collapse of a neutron star to its hybrid star mass twin on the ’’third family” branch, using a 
recent equation of state of dense matter. The available energy corresponding to the mass-energy 
difference between configurations is comparable with energies of the most violent astrophysical 
burst processes. An observational outcome of such a dynamical transition might be fast radio 
bursts, specifically a recent example of a FRB with a double-peak structure in its light curve. 

Subject headings: stars: neutron — radiation mechanisms: non-thermal — gamma-ray bursts: general 
— pulsars: general 


1. Introduction 

Research on the neutron stars’ (NS) equation of 
state (EoS) is currently a very active area. Numer¬ 
ous observations have changed our understanding 
of the cold, dense nuclear matter in compact star 
interiors. Recently, accurate determination of the 
high mass of about 2 M© for PSR J0348-I-0432 
(Antoniadis et al. 2013) and PSR J1614-2230 (De- 
morest et al. 2010) pulsars has provided a power¬ 
ful constraint for the stiffness of the EoS. On the 
other hand, radius measurements are not yet pre¬ 
cise enough to rule out some of the many alter¬ 
native EoS models. Among various estimates, fre¬ 
quency resolved pulse shape analysis for the near¬ 
est millisecond pulsar (Bogdanov 2013) supports 
relatively large radii and offers interesting perspec¬ 
tives for future radius measurements as planned 
with upcoming missions, cf. NICER (Arzouma¬ 
nian et al. 2014). 

Astrophysical models involving NS can ac¬ 
count for the observed powerful gamma ray bursts 
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(GRB, Klebesadel et al. 1973) emission, where 
the typical total energy release is of about lO"*® 
- lO®** erg s“^ (Cavallo & Rees 1978; Paezynski 
1986; Eichler et al. 1989). Other energetic phe¬ 
nomena like fast radio bursts (FRB) - millisecond 
duration radio bursts from seemingly cosmological 
distances (Lorimer et al. 2007), for which a model 
of a NS collapse to a black hole was proposed 
(Falcke & Rezzolla 2014), as well as natal pulsar 
kicks (Podsiadlowski et al. 2005; Lai et al. 2001; 
Kusenko 2005; Berdermann et al. 2005; Stasielak 
et al. 2007) are associated with comparably high 
amounts of energy. 

This Letter considers the energy release during 
a dynamical NS collapse induced by a deconfine¬ 
ment phase transition in the core of a compact star 
(a corequake). Unlike previous works which esti¬ 
mated the energy reservoir for typical (~ 1.4 Mq) 
mass stars (Drago & Tambini 1999; Berezhiani et 
al. 2002; Aguilera et al. 2004; Zdunik et al. 2008), 
we employ a recent EoS model derived in Benic et 
al. (2015) that allows for the formation of a ’’third 
family” of compact stars near the maximum mass. 
For the first time the corequake scenario is consid¬ 
ered in which a high-mass hadronic NS collapses 
into a hybrid compact star disconnected from the 
former by a sequence of unstable configurations 
(for a recent classiheation of hybrid stars, see Al- 
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ford et al. 2013). The point of instability can be 
reached by the hadronic NS in a process of dipole- 
emission spin-down, or accretion-induced spin-up 
by matter from a companion. 

While a direct collapse of a magnetized NS 
to a black hole is one viable mechanism for the 
explanation of FRBs (Falcke & Rezzolla 2014), 
the NS instability and collapse induced by the 
deconfinement phase transition of the type dis¬ 
cussed here would additionally provide an expla¬ 
nation to a double component reported recently 
for FRB121002 by Champion (2015) with a light 
curve peak separation of about 5 ms. In this case, 
the NS would undergo a corequake transition to a 
meta-stable twin hybrid NS configuration, which 
would generate the first peak, before the ultimate 
collapse to a black hole, generating the second 
peak, as in Falcke & Rezzolla (2014) model. 

2. High-mass twin equation of state 

The NS mass twin scenario describes the sit¬ 
uation in which two stars have the same gravi¬ 
tational mass but different radii due to qualita¬ 
tively different internal composition: one of them 
is purely hadronic while the other (the twin) is 
a hybrid star whose core contains quark mat¬ 
ter (Gerlach 1968; Kampfer 1981; Schertler et al. 
2000; Glendenning & Kettner 2000; Dexheimer et 
al. 2015). High mass NS twins (Benic et al. 2015) 
can potentially identify a critical endpoint in the 
QCD phase diagram (Alvarez-Castillo & Blaschke 
2015; Blaschke et al. 2013). At the same time 
they provide an attractive solution to actual prob¬ 
lems of compact star physics discussed in Blaschke 
& Alvarez-Gastillo (2015): the masquerade phe¬ 
nomenon (Alford et al. 2005), the hyperon puzzle 
(Baldo et al. 2003), and the reconfinement prob¬ 
lem (Lastowiecki et al. 2011; Zdunik & Haensel 
2013). An assessment of NS twin identification 
has been recently carried out (Ayriyan et al. 2014; 
Alvarez-Gastillo et al. 2015; Blaschke et al. 2014). 

NS models presented here are based on the DD2 
hadronic EoS with excluded volume correction 
resulting from the quark substructure of nucle¬ 
ons (protons and neutrons) originating from Pauli 
blocking effects at the level of quark substructure 
(Ropke et al. 1986). At suprasaturation densities 
this correction is necessary and has the immediate 
effect of stiffening the EoS. High NS masses are 
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Fig. 1.— Gravitational mass M vs. baryonic mass 
Mb for the hybrid EoS model with vector coupling 
parameter = 0.0. Dynamical NS instability in¬ 
duced by a strong first order phase transition: the 
last stable configuration on the hadronic branch 
(initial state) is connected with the final hybrid 
star configuration at the same baryonic mass with 
a dashed red arrow line. Length of the line rep¬ 
resents the mass difference (energy reservoir) oc¬ 
curring in this transition. The results for 774 =5.0, 
10.0 are qualitatively and quantitatively similar. 

reached already at moderate densities but large 
radii (~ 13 — 15 km, Blaschke et al. 2013). At 
higher densities this EoS undergoes a phase tran¬ 
sition to quark matter described by a NJL EoS 
with multiquark interactions first introduced in 
Benic (2014). The coupling strength parameter 
in the vector channel of the 8 -quark interaction 
( 774 ) determines a sufficient stiffening of the quark 
matter EoS at high densities to preserve the max¬ 
imum observed star mass of 2 Mq in this class of 
models. 

3. Results and discussion 

The EoS of Benic et al. (2015) features the 
possibility of first order phase transitions sub¬ 
stantial enough to destabilize the NS. The ap¬ 
pearance of a small, sufficiently dense quark core 
makes some configurations unstable against radial 


2 




oscillations (see e.g. Shapiro & Teukolsky 1986) 
and is therefore not realized in astrophysical set¬ 
tings. On the mass-radius M{R) diagram for non¬ 
rotating stars these configurations are character¬ 
ized by dM/dXc < 0, where Ac is an EoS parame¬ 
ter labelling the configurations (e.g., their central 
density). Figure 1 traces the subsequent evolution 
- a dynamical collapse - of an unstable compact 
star from the point of the maximum mass Mmax.h 
on the hadronic branch to the corresponding one 
on the hybrid star branch. For simplicity we as¬ 
sume that the baryon mass Mb is conserved dur¬ 
ing the collapse. Figure 2 summarizes the radius 



Fig. 2.— Mass difference AM (upper panel) and 
radius difference Ai? (lower panel) resulting from 
an instability and subsequent NS collapse induced 
by a deconfinement phase transition for a set of 
vector coupling parameters 774 of the high-mass 
twin EoS models. An associated energy release 
AE is indicated on the right side of the upper 
panel. 

change Ai? in the transition, which is between 1 
km and 1.5 km for the coupling constant range 
considered. The energy reservoir AE available in 
the transition, defined as the mass-energy differ¬ 
ence AM(P between the initial and final config¬ 
urations, amounts to roughly 10®^ erg. We note 
that in case of the EoS employed in Zdunik et al. 
(2008) the energy release depends weakly on the 
rotation rate; rotating configurations and the de¬ 


pendence of the energy reservoir on the rotation 
rate for EoS of Benic et al. (2015) will be studied 
in a subsequent article (Bejger et al. 2015). 

In conclusion, the energy reservoir available in 
a high-mass NS instability induced by the decon¬ 
finement phase transition is comparable with val¬ 
ues measured in burst phenomena associated with 
most energetic astrophysical processes. It is clear 
that the deconfinement phase transition in com¬ 
pact stars might play an important role, perhaps 
as an engine for GRBs or as an additional mecha¬ 
nism supporting the explosion of core collapse su¬ 
pernovae, with an accompanied characteristic neu¬ 
trino signals. A magnetized NS that is collapsing 
in a dynamical timescale (~ 1 ms) to another, 
compact configuration may also be attractive in 
the context of FRBs, specifically in view of recent 
observations of FRB121002 (Champion 2015) with 
a double peak light curve structure. Finally, we 
note that the NS EoS is a key input for scenarios 
of cosmic ray generation like supernova explosions 
and NS mergers. Further studies to expose the de¬ 
tails of astrophysical observables in the case of NS 
instabilities induced by dense-matter phase tran¬ 
sitions are work in progress (Bejger et al. 2015). 
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